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TURNOVER  COMPARTMENTALIZATION 

A  Statistical  Tracer  Theory  for  Analysis  of 
Whole -System  Retention  Data 


The  application  of  tracer  methods  have,  especially  in  medicine,  opened 
a  new  class  of  physical  problems,  which  have  seemingly  been  given  little 
consideration.  Probably  this  is  because  the  theoretical  physicists  today  are 
so  much  occupied  by  quantum  mechanics  and  theory  of  relativity  that  there  is 
hardly  any  interest  left  over  for  what  belongs  to  old  fashioned  statistical 
mechanics.  Therefore,  in  the  interpretation  of  tracer  data,  people  have  been 
forced  to  rely  mainly  on  inductive  methods,  i.e.,  the  construction  of  "models" 
that  are  not  derived  from  basic  principles  but  borrow  their  significance  from 
how  well  they  fit  data.  In  an  attempt  to  find  a  more  deductive  approach,  I 
started  some  years  ago  a  construction  of  what  was  intended  to  be  a  general 
physical  tracer  theory,  a  tracer  dynamics  (l).  Although  the  emerging  theory 
was  general  enough,  and  made  possible  a  more  precise  formulation  of  some 
existing  problems,  it  did  not  lead  to  anything  of  much  practical  result.  One 
exception,  however,  can  be  noted:  an  ergodic-like  relation  (P)  which  in  the 
following  shall  be  explained  and  discussed. 

Let  us  first  consider  the  practical  biologic  problem  in  terms  of  a 
specific  example.  In  a  normal  adult  human  there  is  a  daily  uptake  and  a 
balancing  output  of  natural  potassium  ions.  The  result  is  a  fairly  constant 
amount  of  such  ions  in  the  body.  This  naturally  present  substance  (in  the 
present  example,  potassium)  is  in  the  following  referred  to  as  mother 
substance.  A  common  experiment  for  the  study  of  the  potassium  metabolism 
is  to  inject  a  small  amount  of  potassium-42  ions  (the  chloride,  for  instance), 
the  tracer,  and  then  to  measure  how  much  K  remains  in  the  body  at  different 
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times;  frequently  one  also  observes  the  specific  activity  in  blood  (i.e.,  the 
amount  of  tracer  in  a  sample  of  blood  divided  by  the  amount  of  mother  sub¬ 
stance  in  that  sample).  The  problem  then  is:  on  the  basis  of  such  data  what 
can  be  said  about  the  mother  substance? 

One  answer  is  given  by  the  following  equation  (3): 
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where  bg  =  the  total  steady-state  amount  of  mother  substance  in  the  whole 
system  (organism)  S 

=  the  total  amount  of  tracer  in  S  at  the  time  t 
%(t)  =  the  specific  activity  at  the  time  t  in  some  arbitrary 
subsystem  c  S 

The  term  steady  state  means  that  the  system  is  such  that  everywhere 
within  S  the  amount  of  mother  substance  remains  constant  in  time.  The 

physical  significance  of  this  condition  will  be  discussed  in  the  following 
sections. 

The  total  content  of  mother  substance  can  thus  be  obtained  from  this 
type  of  experiment.  However,  there  are  several  critical  remarks  that  could 
be  made  about  equation  1: 

(i)  In  practical  applications  one  would  expect  that  there  are  subsystems 
in  S  with  such  a  low  turnover  that  hardly  any  significant  amount  of 
tracer  ever  appears  there,  although  they  might  contain  a  significant 
amount  of  mother  substance.  These  subsystems  will  not  contribute  to 
the  data  (the  integrals)  and,  therefore,  the  right  hand  member  of 
equation  1  does  not  determine  b°  but  rather  a  smaller  amount, 
commonly  referred  to  as  the  exchangeable  mass  (3).  This  quantity 
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lacks  any  precise  physical  definition  of  its  own,  and  has  to  he 
introduced  as  an  heuristic  concept  (*0  • 

(ii)  It  is  often  difficult,  if  not  even  impossible  (as  seems  to  be  the 
case  for  ^2K+),  to  follow  the  curves  sufficiently  long  to  obtain 
reasonably  good  estimates  of  the  integrals. 

(iii)  Different  functions  bg(t)  can  give  the  same  integral,  which  means 
that  in  equation  1  no  attention  is  paid  to  the  shape  of  the  tracer 
curves.  In  other  words,  the  tracer  curves  must  contain  more  infor¬ 
mation  than  is  revealed  by  equation  1;  i.e.,  the  exchangeable  mass 
cannot  be  the  only  parameter  characterizing  mother  substance  that 
should  be  obtainable  from  tracer  data. 

In  the  following  a  new  approach  to  tracer  analysis  will  be  outlined, 
which  positively  meets  the  objections  (i  -  iii). 

BASIC  PRINCIPLES  OF  THE  APPROACH 


A  rough  formulation  of  the  problem  would  be  that  what  is  needed  is  a  more 
detailed  analysis  of  tracer  curves,  giving  a  more  detailed  description  of  the 
mother  substance,  its  distribution  and  "kinetic  properties"  within  S.  The 
most  commonly  used  approach  today  is  EFCA,  linear  first  order  compartment 
analysis  (5).  The  metabolic  system  S  is  subdivided  into  a  finite  number  of 

so-called  compartments  C^  Cg, . ,CQ,  where  n  usually  is  somewhere  between 

1  to  6.  The  underlying  assumption  is  that  all  tracer  flows  from  Ck(i=l,2, . . .n) 
are  linear  functions  of  the  tracer  content  in  Ck  .  In  other  words,  one  assumes 
first-order  kinetics,  which  occasionally  is  interpreted  such  that  each  Ck 
constitutes  a  "homogeneous"  or  "well- stirred"  subsystem  of  S;  see  further  the 
section  on  the  concepts  of  state  and  homogeneity.  The  compartment  model, 
which  thus  results,  is  represented  by  a  system  of  first-order  linear  differ¬ 
ential  equations  with  constant  coefficients,  having  the  tracer  amounts  in  the 
compartments  as  t -functions.  By  fitting  this  equation  system  to  the  observed 
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functions  bg(t)  and  'a^(t)  one  obtains  estimates  of  the  coefficients,  which 
then  give  the  fluxes  of  mother  substance  (steady  state  assumed)  between  the 
compartments  and  to  the  environment  of  S,  as  well  as  the  amounts  of  mother 
substance  in  the  different  compartments. 

However,  the  concept  of  homogeneity  applied  to  the  compartment  concept 
implies,  as  will  be  discussed  later,  some  basic  difficulties  that  obscure  the 
physical  significance  of  LFCA  (5,6),  and  another  approach  has  been  suggested 
by  Marshall  (7 ) .  Instead  of  considering  the  metabolic  pool  as  subdivided  into 
compartments  defined  in  morphological  and  chemical  terms,  Marshall  considers 
compartmentalization  in  terms  of  "turnover  times."  But  Marshall  works  with  a 
specific  system  (calcium  metabolism),  and  his  method  does  not  seem  to  be 
generally  valid.  In  the  following  I  shall  outline  a  modification  of  Marshall's 
approach  that  seems  to  possess  the  desired  generality.  The  emerging  theory 
has  a  solid  physical  significance,  based  on  earlier  works,  and  the  resulting 
representation  is  applicable  to  a  broad  class  of  systems,  whether  or  not  they 
are  homogeneous  (well-stirred).  The  applicability  is  not  only  in  metabolic 
tracer  work,  but  also  in  biochemistry  (in  vitro  kinetics  of  complicated 
enzymatic  systems)  as  well  as  in  technologic  tracer  methodology. 


THE  ERGODIC  PROPERTY 


I  shall  try  to  limit  the  use  of  abstract  arguments  to  a  minimum.  The 
essential  parts  of  the  mathematical  formalism  have  already  been  worked  out  (8). 
I  shall  use  the  term  particle  meaning  any  kind  of  physical  object  that  retains 
its  identity  throughout  the  system  S •  the  "particle"  may  thus  be  an  ion,  a  dye 
molecule,  a  labeled  fish  and  so  on. 

The  theory  is  based  on  the  following  postulate :  While  in  S  the  particle 
is  at  any  time  in  a  definite  state.  By  this  I  mean  that  at  any  time  a  complete 
description  of  the  particle's  physical  condition  is  in  principle  possible; 


this  does  not  imply  that  such  a  laborious  procedure  has  to  he  carried  out. 

If  the  particle  is  in  some  state  x^  at  the  time  t^,  there  exists  a  definite 
probability  for  this  particle  to  be  in  some  other  state  Xg  at  the  time  tg. 
This  probability  depends  only  on  3^,  Xg  and  tg-t^  and  is  thus  the  same  for 
all  particles  and  all  values  of  it  has  finite  derivatives,  of  all  orders, 
in  the  difference  tg-t^ 

In  a  more  sophisticated  language  the  postulate  implies  that  the  movement 
of  a  particle  through  S  is  a  stationary,  time -homogeneous  Markovian  process, 
but  there  is  nothing  important  in  that  terminology  and  we  can  instead  concen¬ 
trate  our  attention  on  Figure  1.  The  rectangle  denotes  a  closed  S,  and  each 
possible  position  within  it  represents  a  definite  state  of  a  particle;  it  is 
assumed  that  no  matter  what  state  the  particle  is  in  it  can  sooner  or  later 
always  reach  any  other  state  within  S  (the  system  is  said  to  be  irreducible). 


Figure  1 


The  particle  moves  thus  at  random  around  in  S.  Assume  that  at  the  time 
t  =  0  we  observe  a  particle  in  position  (state)  0,  and  then  follow  its  move- 
ments  until  it  finally  reaches  the  state  B  at  t  =  T.  As  illustrated  in 
Figure  1,  the  particle  has  not  passed  through  all  states,  some  states  only 
once,  other  states  twice  and  so  on.  If  now  T  increases,  the  particle  will 
eventually  have  passed  through  all  the  states  (9) >  but  some  states  more  often 
than  others.  Assume  that  T  is  very  large,  and  that  we  at  some  arbitrary  time 
t  >  T  observe  the  particle  in  the  state  x:  the  number  of  trajectories  through 
x  (i.e.,  the  number  of  times  the  particle  has  passed  through  x)  can  then  be 
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used  to  measure  the  probability  of  the  particle  to  be  found  in  that  particular 
state. 

This  has  the  following  important  consequence.  Consider  the  situation 
where,  instead  of  one  single  particle,  S  contains  a  very  large  number  of 
identical  particles,  strolling  around  independently  of  each  other.  It  can 
then  be  shown  (10 )  that  if  the  system  is  left  undisturbed  for  a  sufficiently 
long  time,  it  will  eventually  come  to  a  stationary  state,  or  equilibrium, 
which  means  that  the  density  (or  number)  of  particles  in  each  state  remains 
constant  in  time.  Moreover,  and  this  is  the  very  crucial  thing  indeed,  it  can 
be  shown  that  for  any  state  x  inside  S,  this  equilibrium  density  of  particles 
is  directly  proportional  to  the  number  of  times  our  previous,  lonely  particle 
has  passed  through  that  state. 

What  makes  this  important  is  that  by  observing  a  time  process  (the  move¬ 
ment  of  a  particle)  we  can  estimate  the  stationary  (i.e.,  time  independent) 
distribution  of  particles.  Commonly  this  is  referred  to  as  an  ergodic  property, 
and  the  kinds  of  system  here  considered  (closed  and  irreducible)  are  therefore 
frequently  called  ergodic  systems. 

Let  us  now  consider  an  open  system,  i.e.,  a  system  which  a  particle  will 
eventually  leave  no  matter  what  states  it  passes  through  (Fig.  2).  To  simplify 
the  discussion  also  this  S  is  assumed  to  be  irreducible  (there  is  no  serious 
loss  of  generality  by  this  assumption).  It  is  clear  that  a  particle  entering 
S  will  then  spend  only  a  finite  time  in  S,  which  means  that  there  is  no  justi¬ 
fication  in  assuming  T  to  be  arbitrarily  large.  Consequently,  the  trajectory 
described  by  this  particle  (the  solid  one  in  Figure  2)  can  hardly  be  expected 
to  pass  through  all  possible  states  in  S:  observations  of  the  particle  will 
therefore  not  be  representative  of  the  whole  system.  However,  if  we  repeat 
the  experiment  with  another  particle,  the  obtained  trajectory  will,  with  large 
likelihood,  be  different  (the  dotted  trajectory  in  Figure  2).  Thus,  if  we 
repeat  the  experiment  over  and  over  again,  we  shall  eventually  obtain  a  family 
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of  trajectories  that  covers  the  whole  S.  It  can  then  be  shown  that  the  number 
of  trajectories  passing  through  some  state  x  in  S  has  definite  physical 
significance  (ll): 


Figure  2 


Let  n(x)  be  the  number  of  trajectories  passing  through  x.  Assume  that 
there  is  a  constant  input  q°  of  particles  to  S  (q°  =  number  of  particles 
entering  S  per  unit  of  time).  It  is  then  easily  shown  (12)  that  if  this 
situation  prevails  for  a  sufficiently  long  time,  the  system  will  come  to  a 
steady  state :  the  input  of  particles  becomes  equal  to  the  output  and,  moreover, 
the  number  of  particles  in  each  state  x  remains  constant  in  time.  If  b°(x) 
denotes  the  steady-state  amount  of  particles  in  the  state  x,  the  following 
relation  can  be  proved  (13 ): 

b°(x)  =  q°k  n(x)  Eq.  2 

where  k  is  a  positive  constant  to  be  defined  later.  In  other  words,  the 
steady-state  density  of  particles  in  a  state  x  in  S  is  proportional  to  the 
number  of  tracer  trajectories  through  that  state. 

The  similarity  between  equation  2  and  ergodic  property  that  was  formulated 
for  the  closed  system  is  obvious:  equation  2  may  therefore  be  referred  to  as 
an  ergodic  relation  for  the  open  system,  and  is  identical  with  the  one  previ¬ 
ously  mentioned  as  the  origin  of  equation  1,  which  in  the  following  will  be 
the  main  topic  of  discussion. 

Before  proceeding,  however,  let  us  first  notice  and  appreciate  the 
fundamental  differences  between  open  and  closed  S.  When  dealing  with  a  closed 
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S,  we  were  considering  a  single  particle,  its  trajectory  within  S  during  the 
time  T,  and  we  could  carry  out  a  limiting  process  (i.e.,  allow  T  to  increase 
to  arbitrarily  large  values);  the  number  of  trajectories  through  a  state  x 
will  then  measure  the  probability  of  finding  that  particular  particle  in  that 
particular  state.  Nothing  similar  is  possible  for  an  open  S,  but  we  have  to 
repeat  the  experiment  over  and  over  again,  and  interpret  the  number  n(x)  of 
trajectories  through  x  not  as  the  probability  of  a  single  particle  but  in 
terms  of  a  large  population  of  particles  in  steady  state. 

Moreover,  when  considering  a  closed  S  we  made  use  of  an  equilibrium, 
stationary  distribution.  An  open  S,  however,  can  never  achieve  such  a  state; 
the  ordinary  equilibrium  of  closed  S  is  thus,  for  an  open  S,  replaced  by  the 
steady  state.  That  is,  for  a  closed  S,  the  stationary  distribution  is  the 
natural  reference  system  (time— independent)  in  terms  of  which  we  can  interpret 
the  observed  time-process,  and  for  an  open  S  the  steady-state  distribution  is 
the  corresponding  reference  system.  This  is  the  reason  why,  when  considering 
open  metabolic  systems,  the  steady-state  assumption  is  absolutely  necessary! 

(iM- 


THE  APPLICATION  TO  TRACER  THEORY 

The  application  is  straightforward.  We  can,  in  equation  2,  identify  q° 
with  the  input  (daily  uptake)  of  mother  substance,  b°(x)  with  the  amount  of 
mother  substance  in  the  state  x  when  S  is  in  steady  state,  and  n(x)  with  the 
number  of  times  tracer  particles  (no  matter  which  ones)  have  passed  through  x. 
The  remaining  problem  is  the  factor  k  and  it  is  solved  as  follows; 

n(x)  can  be  interpreted  as  being  directly  proportional  to  the  time  of 
sojourn  9(x)  in  state  x,  i.e.,  the  total  time  a  particle  is  "expected"  to 
spend  in  x.  In  other  words 


n(x)  =  k'  9(x) 
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where  kf  is  another  constant  independent  of  x.  If  we  consider  an  arbitrary 
subsystem  (domain)  in  S,  the  sum  over  all  states  x  in  is 


=  ^  k'9(x)  =  k' 


where  6  is  the  time  of  sojourn  in  S  ,  i.e.,  the  total  time  a  particle  is 
x>  u 

expected  to  spend  in  S  .  Summing  over  all  states  in  S, 


ng  =  Y  k’e(x)  =  k 


which  thus  gives  the  total  time  a  particle  is  expected  to  spend  in  S,  i.e., 
between  entrance  and  exit  (15)*  Finally  it  can  be  shown  (l6)  that  k.kf  =  1 


and  that 


(x)  =  doT  I  b(x't} 


where  b(x,t)  is  the  amount  of  tracer  in  x  at  the  time  t,  and  bg(o)  is  the  total 
amount  tracer  supplied  to  S  at  the  start  of  experiment  (t  =  0).  It  is  thus 
assumed  that  there  is  a  single  instantaneous  tracer  supply,  but  this  assumption 
is  not  necessary  and  is  introduced  merely  to  simplify  the  formalism. 

Using  these  results  we  obtain  from  equation  2  that 


=  =  TJO 

o 


J 


where  h°  is  the  steady  state  amount  of  mother  substance  in  S  and  b  (t)  is 
*0  X)  v 

the  corresponding  amount  of  tracer  at  the  time  t.  Writing  the  subscript  t>  =  S 
we  get  the  same  quantities  for  the  whole  system  (i.e.  S  =  S)  and  if  we  define 


ap(t)  = 


equation  1  is  immediately  obtained. 
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My  motivation  for  using  this  particular  way  of  deriving  equation  1  is 
that  it  provides  us  with  a  tool,  by  which  we  can  easily  visualize  the  physical 
significance  of*  the  basic  equation 

,  o  o^ 

bU  =  *  %  Eq.  k 

namely  that  the  steady-state  amount  of  mother  substance  in  some  arbitrary 
subsystem  is  proportional  to  the  density  of  tracer  trajectories  through 
that  domain.  This  gives  us  an  intuitive  understanding  of  why  and  how  tracer 
tells  us  something  about  mother  substance.  And  it  is  on  the  basis  of  this 
interpretation  of  equation  k  we  shall  be  able  in  the  next  section  to  make  use 
of  the  previously  mentioned  idea  induced  by  Marshall. 

Before  concluding  this  section  there  is  another  quantity  that  has  to  be 
introduced.  Let  us  consider  a  particle  in  some  state  x.  As  long  as  the 
particle  remains  there  it  has,  according  to  the  postulate,  a  definite  proba¬ 
bility  to  leave  the  system  S  within  some  fixed  time  interval.  Denote  this 
probability  by  X£(x)and  consider  the  expression 

XEv  =  7Z  I  b°(x )  Vx) 

b  s 
v 

Hence,  is  the  mean  exit  probability  in  S  ,  calculated  with  respect  to  the 
steady-state  distribution  in  that  subsystem.  By  putting  t>  =  S  we  obtain  the 
mean  exit  probability  in  the  whole  system  with  respect  to  the  steady-state 
distribution  of  mother  substance  in  S.  It. holds  that  (17) 

^ES  =  Veg  Eq.  5 

Hence,  as 

00 

®S  =  bjoj  J  bS(t)  dt  E**  6 

S  o 

the  area  under  the  'whole-body  retention  curve  can  be  given  a  precise  micro- 

physical  interpretation:  it  determines  the  mean  exit  probability  A°  with 

ES 

respect  to  the  steady-state  distribution. 
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To  me  this  is  an  important  result,  because  it  shows  that  data  from 
systems  of  enormous  microphysical  complexity  can,  occasionally,  he  given  a 
precise  microphysical  interpretation.  In  some  way  is  analogous  to  the 
temperature  concept  in  physics:  temperature  can  he  interpreted  as  a  parameter 
characterizing  the  mean  energy  over  the  population  of  molecules  (i.e.,  temper¬ 
ature  is  given  a  microphysical  significance).  But,  for  this  interpretation  to 
he  valid  the  physical  system  must  he  in  thermal  equilibrium  (the  Bolzman 
distribution,  for  instance),  and  similarly  for  the  interpretation  given  here 
of  retention  curves,  S  must  he  in  steady  state.  So  here  again  we  have  an 
illustration  of  how  the  steady-state  distribution  plays  a  role  for  the  open 
system  similar  to  that  of  the  equilibrium  or  stationary  distribution  in  the 
theory  of  closed  systems. 


TURNOVER  COMPARTMENTS 

There  is  one  major  conclusion  that  can  he  made  from  the  preceding  analysis: 
it  is  the  complete  tracer  process,  and  not  separate  parts  of  it  that  tells  us 
something  about  mother  substance;  this  is  why  integrals  appear  in  the  equations. 
Phrased  differently,  it  can  be  said  that  one  has  to  follow  the  trajectory  of  a 
tracer  particle  from  its  beginning  to  its  end  (i.e.,  when  the  particle  leaves 
S).  However ,  one  single  trajectory  is  a  random  event  and  gives  a  poor  or, 
rather,  uncertain  information  about  the  system.  It  is,  therefore,  necessary 
to  carry  out  a  large  number  of  statistically  independent  "experiments”,  a 
condition  that  actually  is  rather  well  satisfied  in  a  tracer  experiment :  each 
tracer  particle  represents  one  independent  trial.  A  tracer  experiment  might 
be  looked  upon  as  generating  a  large  population  of  trajectories  through  S,  and 
the  observed  tracer  variables  may  be  considered  as  averages  over  this  popula¬ 
tion,*  the  magnitude  of  AvogadroTs  number  implies  that  the  statistical  error  is 
negligible . 
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It  is  an  obvious  statement  that,  among  the  population  of  trajectories 
generated  by  a  tracer  experiment,  some  are  fairly  short  whereas  others  are  of 
a  considerable  length.  That  is,  some  tracer  particles  leave  S  shortly  after 
they  entered,  whereas  other  particles  will  stay  in  the  system  for  quite  a 
while.  Hence,  we  can  subdivide  the  trajectory  population  into  subpopulations 
with  respect  to  the  length  of  the  different  trajectories. 

Let  us  consider  the  total  t-interval  (0,T)  during  which  the  experiment  is 
carried  out,  and  divide  it  into  intervals 


0  <  tl  <  t2 


<  t  ,  <  t 
n-1  n 


=  T  . 


As  before,  bg(t)  denotes  the  amount  of  tracer  present  in  S  at  the  time  t;  i.e., 

the  whole-body  retention  curve.  Then  bg(o)  -  bg (t^)  is  the  number  of  tracer 

particles  leaving  S  in  the  first  t-interval  and,  generally,  b  (t.  )  -  b  (t.) 

S  i— 1  S  i 

is  the  number  of  particles  leaving  S  in  the  ith  interval.  This  corresponds 
to  a  subdivision  of  the  total  population  of  trajectories,  as  discussed  in  the 
preceding  paragraph.  But,  as  the  density  of  trajectories  in  each  such  subpopu¬ 
lation  corresponds  to  a  steady-state  density  of  mother  substance,  the  procedure 
simply  implies  a  subdivision  of  the  total  pool  of  mother  substance  with  respect 
to  when  the  particles  leave  S.  Consequently,  we  have  a  division  of  the  popula¬ 
tion  of  mother  substance  particles  into  n  mutually  exclusive  turnover 
compartments  K^,  K^,  K^,  each  of  which  has  a  definite  steady-state  amount 

of  mother  substance  and  a  "kinetics"  (i.e.,  a  -parameter)  of  its  own.  The 
following  will  show  how  these  quantities  can  be  estimated  from  b  (t)-curves. 

Since  in  the  present  context,  the  main  purpose  of  the  integral  of  a^(t) 
in  equation  1,  is  to  determine  q°  this  integral  will  from  now  on  be  omitted: 
the  amount  of  mother  substance  is  always  directly  proportional  to  q°,  and  we 

shall  thus  here  satisfy  ourselves  with  the  relative  sizes  of  KL,  . .  K  If 

- —  Tm  n 

the  a(t) -integral  can  be  obtained,  the  absolute  sizes  are  easily  calculated 
from  the  quantities  obtained  by  the  methods  described  in  the  following. 
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Let  us  start  by  defining  the  function 


P.(t)  =  ■  y 

1  bg(0)  -  bg(t^)  i  =  1,2,..., n 


Eq.  T 


Consequently,  P.(0)  =  1,  Pjq)  =  0,  and  the  function  itself  may  be  interpreted 
as  the  probability  that  a  particle  will  be  in  S  at  the  time  t,  given  that  it 
belongs  to  any  of  the  compartments  K^,  Kg,  ...  K^.  That  is,  EL (t)  is  a  so- 
called  conditional  probability,  where  the  condition  is  that  the  particle  is 


known  to  leave  S  somewhere  during  the  time  t.  (we  start  at  t  =  0). 


Eq.  8 


is  the  conditional  time  of  sojourn.  To  get  the  "absolute  time  of  sojourn 

0  0*  must  be  multiplied  by  the  probability  that  a  particle  entering  S 

o,i'  o,i 

will  belong  to  the  set  of  particles  corresponding  to  the  joint  compartments 
K^,  Kg  ...  K1>  This  probability  is  given  by  [bg(0)  -  bg(ti)]/bg(0)  and. 


consequently 


e  . 

0,1 


*i 

I 


bs(t)  -  bg(t . ) 

- -  dt 

bg(0) 


Eq.  9 


Finally,  if  9^  denotes  the  time  of  sojourn  of  K^  alone. 


e.  =  e  .  -  e  .  , 

1  0,1  0,1-1 


Eq.  10 


Here  0.  has  the  same  significance  as  0^  in  equation  k .  That  is 

h°  =  q°.0.  Eq.  11 

l  i 

is  the  steady-state  amount  of  mother  substance  in  compartment  in  other 
words  b°  is  the  ’’size”  of  K±,  and  0i  determines  what  we  then  call  the  relative 

size  of  K^# 
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Compartment  can  be  meaningfully  subdivided  into  two  mutually  exclusive 
subcompartments : 


=  the  initial  -compartment 
Ki^  =  the  final  -compartment 

where  corresponds  to  those  particles  in  that  have  spent  a  time  less 
than  an(^  "those  particles  in  that  have  spent  a  time  greater 

than  t_^^.  Only  contributes  to  the  output  of  S;  i.e.,  the  exit  probability 
of  particles  in  is  zero. 

From  equation  9  we  understand  that,  because  times  of  sojourn  of  mutually 

exclusive  compartments  are  additive,  the  integral  between  t.  and  t  gives 

l-l  i 

the  time  of  sojourn  in  K^,  that  is 


r  -  W 

ei2  =  J  — - S-L- dt 

t±_i  bfO) 


Eq.  12 


and  obviously, 


6.n  =  0.  -  0.o 
ll  l  i2 


Eq.  13 


which  is  the  time  of  sojourn  in  K.^. 

eil  and  e-j_2  meas^re  the  relative  sizes  of  the  corresponding  compartments. 

What  remains  is  to  find  a  measure  of  the  kinetics  of  K. .  However,  as  K  does 

i  9  il 

not  contribute  anything  to  the  output  of  S,  we  are  only  interested  in  K.^, 

Consequently,  we  are  interested  in  obtaining  the  mean  exit  probability  with 

respect  to  the  steady-state  distribution  in  K^.  Thus,  given  that  the  particle 

is  in  K  its  probability  of  still  being  in  S  at  the  time  t  (t  <  t  <  t  ) 

±  i-1  —  —  i' 

is  given  by 


pi2(t)  = 


bs(t)  -  W 

Wi)  -  w\  =  ° 


Eq.  l4a 
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and 


Eq.  14b 


The  corresponding  mean-exit  probability  is  given,  finally,  by 


^°i  =  V4  • 


Eq.  15 


It  is  important  to  recognize  that  0^  does  not  determine  the  size  of  Ki2; 
it  gives  only  the  time  of  sojourn  in  under  the  condition  that  the  particle 
arrives  at  K  according  to  the  distribution  over  the  states  in  K.  at  the  time 

i2  1 

t  :  under  the  present  conditions  only  a  fraction  of  the  total  input  q 
i~l7 

contributes  to  K±.  On  the  other  hand,  \°^(q°0^2)  gives  the  contribution  of  K± 
to  the  total  output  from  S. 


THE  TURNOVER  DIAGRAM 


It  is  essential  that  the  results  of  the  calculations  (the  data  processing) 
can  be  exhibited  in  a  neat  form,  so  that  the  qualitative  interpretation  is 
simplified.  Each  turnover  compartment  (i  =  l,2,...,n)  is  characterized  by 
four  parameters:  0  ,  ©i2,  and  the  time  interval  [t^,  t±] .  Figure  3 

shows  one  possible  type  of  diagram:  the  area  of  each  "box"  gives  the  corre¬ 
sponding  6-value  (the  relative  size),  together  with  the  actual  number  printed 
inside  the  box.  The  length  of  the  arrows  is  proportional  to  XEi  ei2>  1,e*' 
proportional  to  the  relative  contribution  of  each  compartment  to  the  total 
output  of  mother  substance  from  S. 

Figure  3 
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An  advantage  of  the  suggested  diagram  is  that  it  clearly  reveals  the 
fundamental  character  of  the  concept  turnover  compartment:  There  is  no 
a  priori  assumed  relation  ‘between  the  turnover  compartments  and  the  physico¬ 
chemical  structure  of  the  metabolic  system.  If  any  such  relation  exists  it 
has  to  be  shown  by  other  means,  independent  of  the  tracer  experiment.  That 
is,  the  turnover  compartment  is  defined  independent  of  "location”  and  physico¬ 
chemical  "functions",  and  from  that  viewpoint  the  concept  becomes  rather 
abstract.  It  is  also  clear  that  it  makes  no  sense  to  speak  of  exchange  of 
matter  between  turnover  compartments. 


THE  CONCEIT  OF  "STATE" 


A  theory  of  the  type  considered  here  may  be  looked  upon  as  language 
consisting  of  a  set  of  basic  concepts  (words),  and  a  postulate  (grammar) 
ascribing  some  general  properties  to  them.  On  this  basis  all  other  concepts 
and  relations  (equations)  are  derived^  the  basic  concepts  are  thus  not  defined 
but  considered  "obvious".  The  basic  concepts  here  used  are,  particle,  state, 
and  probability  (l8).  Among  these  I  feel  that  the  concept  of  state  needs 
some  further  considerations. 

Let  us  start  with  an  impossible  type  of  system:  S  consists  of  one  single 
state.  According  to  the  postulate  all  particles  in  S  have  the  same  constant 
probability  \  of  leaving  S.  Hence,  if  we  start  at  the  time  t  =  0,  the 
probability  Pg(t)  that  a  particle  will  still  be  in  S  at  some  time  t  >  0  is 
given  by  (19) 


ps(t) 


Eq.  14 


Hence,  the  expected  number  of  particles  remaining  in  S  decreases  according  to 
a  simple  "exponential  law". 
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Next,  consider  a  system,  which  might  he  referred  to  as  a  well - stirred 
system  (Figure  4):  a  small  "open"  vessel  equipped  with  some  stirring  device. 
The  system  is  in  steady  state:  there  is  a  constant  inflow  of  plain  water,  an 
equal  outflow,  and  the  amount  of  water  in  the  vessel  remains  constant.  If 
now  a  small  amount  of  tracer  particles  is  added  to  the  vessel  at  the  time 
t  =  0,  and  the  tracer  content  of  S  is  observed  as  a  function  of  time,  it  is  a 
well  established  experience  that  these  data  will  satisfy  equation  l4.  Hence, 
the  well-stirred  S  behaves  as  if  it  consisted  of  one  single  state,  which, 
however,  is  in  bad  agreement  with  the  intuitive  understanding  of  the  state 
concept. 


Figure  4 


Thus,  if  at  some  instant  we  observe  two  particles  simultaneously,  one  at  x^ 
and  the  other  at  x^,  we  should  hardly  consider  them  as  being  in  the  same  state. 
Moreover,  it  seems  rather  obvious  that,  at  the  moment  of  observation,  the 
particle  in  x^  has  a  higher  probability  of  leaving  S  than  has  the  particle  in 
x .  Consequently,  on  the  basis  of  our  intuitive  understanding  of  the  state 
concept,  we  cannot  accept  this  S  as  consisting  of  one  single  state;  the  other 
possibility,  namely  that  all  states  in  S  are  equal  with  respect  to  the  exit 
probability,  must  be  rejected  as  well. 

As  a  matter  of  fact  we  are  here  faced  with  a  kind  of  a  time -scale  problem. 

Let  X  (x,At)  stand  for  the  probability  of  a  particle  in  state  x  (a  position  x 

E 

in  the  vessel  of  Figure  4)  to  leave  the  vessel  within  the  time  interval  At; 

when  X  (x,At)  is  the  same  for  all  x  then  the  process  becomes  a  simple  exponen- 

tial  e“^EL  where  \  =  XE(x)  is  such  that  XE(x,At)  =  XgAt  +  o(At).  However, 

for  X  (x,At)  to  be  independent  of  x  depends  on  the  scale  in  which  At  is 
E 

measured.  For  instance,  if  At  stands  for  a  certain  number  of  milliseconds, 


IT 


\ (x, At)  will  be  different  for  different  x,  whereas  it  becomes  almost  indepen¬ 
dent  of  x  if  At  is  measured  in  hours.  In  other  words,  if  we  ask  for  the 
probability  of  a  particle  to  leave  S  within  10  hours,  say,  we  shall  obtain 
the  same  answer  for  x^  as  for  x^,  whereas  this  is  not  so  when,  in  the  question, 
hours  is  replaced  by  milliseconds  (it  is  here  assumed  that  the  overall  tracer 
process  is  so  fast  that  after  one  hour  almost  all  tracer  particles  have  left 
S).  That  this  has  corresponding  phenomenologic  consequences  has  been  shown 
elsewhere  (20 ) :  when  the  time  scale  of  observation  is  increased  the  system 
behaves  as  if  the  different  states  becomes  lumped  together,  i.e.,  S  behaves 
"simpler"  in  large  time  scale  than  in  a  small  one.  The  phenomenon  may  be 
referred  to  as  a  time  scale  lumping. 

The  commonly  employed  way  out  of  the  dilemma  is  to  state  that  the  distri¬ 
bution  of  particles  over  the  different  states  in  S  is  time  independent.  Such 
a  time  independency  of  distribution  leads  to  the  desired  system  behavior  (21). 
This  is  why  in  classic  chemical  kinetics  the  different  chemical  compounds  in 
the  reaction  can  be  considered  as  "states":  the  necessary  condition  is  that 
the  system  is  in  thermal  equilibrium,  which  means  that  the  energy  distribution 
over  the  population  of  molecules  is  time -independent  (the  Bolzman  distribution). 
In  these  concepts,  then,  one  could  consider  the  concept  of  "efficient  stirring" 
as  a  generalization  of  the  condition  in  chemical  kinetics.  However,  the  reason 
that  the  concept  of  time -invariant  state  distribution  works  in  chemistry  is 
that  the  condition  of  thermal  equilibrium  seems  to  be  well  satisfied  also  for 
rather  fast  chemical  processes  (22).  In  other  words,  within  the  range  of  time 
units  in  which  chemical  reactions  are  usually  studied,  the  underlying  concepts 
are  invariably  applicable:  what  is  thermal  equilibrium  for  a  fast  reaction  is 
exactly  the  same  thing  for  a  slow  one. 

In  the  present  investigation  we  are,  however,  faced  with  a  different 
class  of  systems.  First,  the  spectrum  of  time  scales  and  experimental  preci¬ 
sion  (absolute  and  relative)  in  which  tracer  studies  are  carried  out  has  a 
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considerable  width.  Second,  the  discussion  of  Figure  k  shows  that  the  well- 
stirred  system  is  highly  heterogeneous  when  considered  in  terms  of  the  two 
basic  concepts,  state,  and  probability.  Due  to  time  scale  lumping  the  phenom- 
enologic  behavior  of  such  a  system  is  likely  to  depend  on  the  method  of 
observation  in  terms  of  time  scale  and  precision  and,  consequently,  what  in 
one  tracer  study  appears  to  be  a  homogeneous  S  (all  states  are  equal)  may  not 
necessarily  be  so  in  another  tracer  experiment.  Or  phrased  differently,  when 
we  consider  the  concept  of  homogeneity  in  terms  of  time -invariance  in  the 
distribution  of  particles  over  the  states,  in  S,  what  under  certain  circum¬ 
stances  appears  to  be  a  time  invariant  distribution  can  under  other  conditions 
show  significant  variation. 

The  concept  of  homogeneity  (or  efficient  stirring),  which  is  further 
discussed  in  the  next  section,  has  played  a  fundamental  role  in  current  tracer 
theories  (23),  where  it  is  usually  defined  in  phenomenologic  terms  (the 
"exponential  law").  Thus,  the  concept  of  the  homogeneous  compartment  is  the 
cornerstone  in  the  LFCA-method  previously  mentioned,  but  as  the  preceding 
paragraph  indicates,  the  approach  is  likely  to  lead  to  a  rather  uncontrolled 
situation:  although  different  investigators  might  use  the  same  terminology 
(compartment,  efficient  stirring,  and  so  on)  the  microphysical  significances 
of  the  words  can  be  quite  different.  The  situation  is  similar  to  describing 
a  moving  body  in  terms  of  a  coordinate  system  that  changes  unpredictably  with 
time.  That  a  similar  approach  works  in  classic  chemical  kinetics  is  due  to 
the  particular  circumstances  there. 

From  this  vie-wpoint  it  is  therefore  important  to  recognize  the  signifi¬ 
cance  of  the  state  concept  introduced  here:  it  is  an  invariant  microphysical 
entity,  which  has  the  same  meaning  for  all  systems  (experiments)  to  which  the 
theory  applies.'  Thus,  whereas  the  number  of  exponential  terms  necessary  to 
describe  a  tracer  curve  is  a  crucial  point  for  many  LFCA -methods  the  present 
theory  puts  no  significance  in  this  system  property.  The  well-stirred  system 
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is  represented  by  the  same  kind  of  turnover  diagram  as  is  the  highly 
"heterogeneous”  system. 


THE  CONCEPT  OF  HOMOGENEITY 

So  far  hg(t)  has  been  considered  both  as  the  actual  function  observed  by 
the  experimenter  and  as  the  result  of  the  observation.  In  the  following  it 
will  be  necessary  to  discriminate  between  those  two  quantities.  Thus,  from 
now  on  bg(t)  denotes  the  physical  f unction  (or  observed  function),  which  is 

defined  by  S  alone  independently  of  experimental  procedures,  i.e.,  ba(t)  is 

b 

defined  only  by  the  set  of  states,  corresponding  to  S,  and  a  transition  proba¬ 
bility  function  on  this  set.  The  experimenter  observes  b  (t)  and  records  the 
observations  in  a  table  or  diagram  and,  by  some  fitting  procedure,  he  ends  up 
with  a  mathematical  form  bg(t),  which  conveniently  describes  the  result  of  the 
complete  experimental  procedure  (including  the  data  processing).  I  shall  call 
bg (t )  the  resulting  function. 

We  may  thus  consider  a  procedure  C  operating  on  the  physical  function 
bs(t)  giving  the  result  (or  image)  b*(t) : 

C[bg(t)]  =  b*(t)  Eq.  15 

There  is  no  reason  to  assume  that  the  mathematical  form  of  hg(t)  is  preserved 

by  the  mapping  C.  Thus,  if  we  still  assume  that  S  contains  a  finite  number  of 

states  bg(t)  will  consist  of  a  sum  of  a  finite  number  of  exponentials,  i.e., 

-Ctf  •  t 

terms  of  the  form  p^e  1  ,  where  01  is  positive  definite  and  t -independent, 
whereas  p±  may  contain  sine  and  cosine  functions  of  t.  It  can  then  be  shown 

■X*  -$f 

(22)  that  bg(t),  written  as  a  sum  of  p^e  ^  (where  p^  is  t-independent),  will 
in  all  likelihood  contain  a  less  number  of  exponentials  than  bg(t);  so  in 
these  concepts  there  is  no  simple  relation  between  bQ(t)  and  b*(t).  It  is 
thus  possible  that  the  resulting  function  contains  only  one  single  exponential. 
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while  the  observed  function  contains  a  fairly  large  number  of  such  terms.  The 
only  condition  under  which  we  can  make  a  rather  safe  statement  about  bg(t), 
on  the  basis  of  the  resulting  function  b*(t),  is  when  bg(t)  consists  of  one 
single  exponential. 

DEFINITION:  S  is  said  to  be  homogeneous  if  and  only  if  bg(t)  consists  of  one 
single  exponential.  S  is  said  to  be  approximately  homogeneous  when  it  is  non- 
homogeneous  but  bg(t)  consists  of  one  single  exponential. 

For  illustrative  purposes  I  shall  discuss  the  simple  two-state  system 
shown  in  Figure  5.  q?  is  the  number  of  mother-substance  particles  that  enter 
S  at  state  i  (i  =  1,2)  per  unit  of  time,  and  {>  }  are  the  "rate  constants" 

or  "probability  factors".  The  instantaneous  tracer  supply  at  t  =  0  is 
throughout  assumed  to  be  equivalent  (24):  the  amounts  of  tracer  added  to  1 
and  2  are  proportional  to  q°  and  q^  respectively.  It  is  only  under  this 
condition  that  tracer  represents  mother  substance,  and  the  ergodic  relation 
holds. 

It  can  easily  be  shown  (25)  that  S  becomes  homogeneous  in  the  present 

sense  when  =  X_0.  Efficient  stirring,  on  the  other  hand,  gives  rise  to 
El  E2 

an  approximately  homogeneous  S,  as  has  been  discussed  in  the  preceding  section. 

Commonly  the  concept  of  efficient  stirring  is  interpreted  as  X^9  »  ^1^2' 

However,  this  condition  does  not  approximate  the  homogeneous  S  defined  by 

X  =  \  ,  and  the  question  becomes  whether  the  behavior  of  the  efficiently 

El  E27 

stirred  S  is  due  to  its  approximation  to  an  homogeneity  condition  and,  if  so, 
what  kind  of  condition  is  it? 

First  of  all,  the  behavior  of  the  efficiently  stirred  S  cannot  be  a 

consequence  of  C  alone,  because  we  are  here  considering  a  property  of  S  itself. 

Second,  as  I  have  pointed  out  in  the  previous  section,  attempts  to  consider 

homogeneity  in  terms  of  time  independency  of  the  tracer  distribution  in  S 

lead  to  time-scale  problems.  We  are  thus  looking  for  a  quantitative  condition 

(other  than  A.—..  =  A.TPO),  which  will  be  insensitive  to  time-scale  variations  (or 
'  El  E2 


21 


the  experimental  procedure  C  in  general)  hut  will  result  in  a  homogeneous  S 
that  efficient  stirring  does  approximate.  The  problem  is  solved  by  the 
following  theorem. 


Figure  5 


1 

"V  2 

THEOREM:  S  is  homogeneous  if 

(i)  S  contains  a  finite  number  of  states 

(ii)  b(x1,0)  =  b(xg,0)  for  all  x^Xg  e  S 

(iii)  h(x^,0)  =  b(x2,0)  for  all  x^Xg  e  S 

That  means,  in  terms  of  Figure  5,  that  q.°  =  and,  if  b^t)  is  the 
amount  of  tracer  in  state  i  (i  =  1,2)  at  the  time  t,  it  also  is  true  that 

VO)  =  V0)-  The  Proof  of  the  theorem  will  show  what  restrictions  on  the 

matrix  [X  }  are  implied  by  the  conditions  (ii)  and  (iii). 


PROOF: 

b(t)  = 


V6) 


b  (t) 
■  nv 


is  the  column  vector  of  tracer  contents  in  the  states 


xq . xn  contains  n  states  altogether.)  For  a  unit  instantaneous  input 

of  tracer  at  t  =  0  (equivalent  input  assumed)  one  has  (26) 

b(t)  =  A  b(t)  [i] 

where  A  =  is  the  matrix  of  transition  probability  factors.  Then 


(ii)  *  b(0)  =  ± 


,  and  (ii)  =>  b(o) 
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Hence,  according  to  [I] 
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Such  a  matrix  A  has  the  eigenvector 


-n.  That  is 


and  the  corresponding  eigenvalue 


A  Vir 


The  solution  of  [I]  can  be  written  as  (27) 

b(t)  =  £  c.  X.  e 


where  X.  is  the  eigenvector  of  A  corresponding  to  the  eigenvalue  a^.  It  is 
assumed  that  A  of  order  n  X  n,  has  n  distinct  eigenvalues  (28).  To  obtain 
the  coefficients  {CL}^=^  we  take  the  first  derivative  of  b(t)  and  calculate 
its  value  at  t  =  0: 


b(0)  =  UI  a.C.  f  x2i 


anC_+a^Xn^C0+ 

....  + 

a  X. 

11  2  12  2 

n  . 

a  C  +a 

....  ~f* 

a  X, 

11  2  22  2 

n  1 

a.Ch+aXX  o0o+ 

....  + 

a  x 

11  2  n2  2 

n  i 

where  we  have  identified  X^  with  the  known  eigenvector 
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The  matrix 


1  a2X12 
1  a2X22 


1  a2Xn2 


a  Xn 
n  In 

a  Xn 
n  2n 


a  X 
n  nn 


is  nonsingular  (29)- 
is  a  solution  to  [V]. 
solution  (29).  Hence 
form 


Moreover,  it  is  easily  seen  that  Cn=  — ,  C  -CL-. .  .=C  =  0 

1  1  2  3  n 

But  as  the  matrix  is  nonsingular,  this  is  the  only 
,  from  [IV]  we  obtain  that  the  solution  of  [I]  has  the 


a^  =  -n 


which  proves  the  theorem. 

Consequently,  the  theorem  states  that  if  a  finite  tracer  system  has  a 
transition  matrix  such  that  the  sum  of  the  elements  in  each  row  gives  the  same 
value  for  all  rows,  then  the  system  will  exhibit  a  single  exponential  tracer 
curve,  provided  the  initial  distribution  of  tracer  is  rectangular  over  the 
states  in  the  system.  In  terms  of  Figure  5*  this  type  of  homogeneous  S  may  be 
illustrated  by  q.°  =  =  1,  X^  =  3,  X£2  =  1,  Xgl  =  1,  \2  =  2.  Certainly, 

this  example  presents  us  with  a  microphysically  heterogeneous  system:  the 
two  states  are  quite  different  in  terms  of  X  and  the  efficient  stirring 
condition  (i.e.,  ^22*^21  »  ^ar  ^rom  satisfied.  I  shall  refer  to 

this  new  class  of  homogeneous  S  as  weak  homogeneity,  as  distinguished  from  the 
strong  homogeneity  defined  by  the  condition  X^  -  (i,j  =  1,2,..., n).  In 
contrast  to  weak  homogeneity,  strong  homogeneity  implies  that  b  (t)  is  a 
single  exponential  no  matter  how  the  tracer  is  added  to  S;  as  a  matter  of  fact, 
strong  homogeneity  implies  that  S  behaves  as  if  it  contained  one  single  state 
(the  different  states  appear  "lumped"  together). 
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Now,  which  class  of  homogeneity  does  the  well-stirred  S  illustrated  by 
Figure  k  approximate?  We  have  already  seen  that  it  does  not  approximate  the 
strong  homogeneity,  and  a  moment’s  contemplation  will  show  that  it  approximates 
the  type  of  weak  homogeneity  defined  by  the  preceding  theorem.  The  fact  that 
its  behavior  is  independent  of  how  the  tracer  is  added  is,  in  this  particular 
case,  irrelevant. 

The  existence  of  weak  homogeneity  is  fundamental  for  our  understanding  of 
the  behavior  of  many  tracer  systems,  and  of  their  striking  simplicity.  Charac¬ 
teristic  of  the  weak  homogeneity  is  that  the  behavior  is  partly  due  to  how  the 
tracer  is  added  to  S,  a  fact  that  again  demonstrates  how  irrelevant  it  is  to 
ascribe  significance  to  the  number  of  exponentials  into  which  bgft)  can  be 
resolved. 


CONCLUDING  REMARKS 

Let  us  now  return  to  equation  1  and  the  criticism  of  it.  Objection  (i) 
is  obviously  met  by  the  approach  I  have  outlined  in  the  preceding  sections. 

The  term  11  exchangeable  mass”  thus  appears  somewhat  inadequate,  at  least  if  it 
is  not  accompanied  by  a  specification  of  the  time  interval  under  consideration. 
Hence,  instead  of  speaking  of  exchangeable  mass  we  now  consider  the  sizes  of 
certain  turnover  compartments,  i.e.,  the  amount  of  mother  substance  that  under 
steady  state  conditions  will  stay  in  S  a  specified  period  of  time. 

A  striking  feature  of  the  theory  is  the  way  in  which  it  meets  objection 
(ii).  Thus,  no  matter  how  short  the  time  period  T  of  observation  is  we  can 
always  construct  a  turnover  diagram.  But,  on  the  other  hand,  the  diagram  ends 
promptly  at  T,  which  means  that  the  theory  does  not  allow  for  any  conclusions 
to  be  made  beyond  the  time  of  observation.  In  other  words,  the  turnover 
diagram  represents  actual  data  and  nothing  else.  This  is  somewhat  different 
from,  say,  LFCA,  which  is  not  a  pure  kinetic  analysis  of  data  in  the  present 
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sense:  opposite  to  the  present  theory,  LFCA  takes  into  account  also  a  priori 
knowledge  or  hypotheses  about  S,  and  makes  possible  a  more  extensive  data 
analysis.  In  my  own  opinion  this  is  a  disadvantage  of  LFCA,  in  that  this 
approach  contains  pitfalls  and  does  not  clearly  discriminate  between  the 
observed  kinetics  of  S  and  specific  hypotheses;  that  is,  there  is  in  LFCA  an 
inbuilt  possibility  to  prove  what  one  wants  to  prove. 

Concerning  objection  (iii)  it  is  clear  that  a  turnover  diagram  can  be 
made  as  detailed  as  one  wishes:  there  is  no  theoretical  limit  as  regards  the 
fineness  of  the  subdivision  of  the  interval  0  to  T.  On  the  other  hand,  there 
is  no  sense  in  making  the  representation  more  detailed  than  is  motivated  by 
data.  Here  also  the  a  priori  knowledge  about  S  becomes  significant:  there  is 
no  reason  for  making  the  turnover  compartmentalization  more  detailed  than  the 
number  of  functions  and  properties  of  S  one  wants  to  study. 

In  this  context  I  should  like  to  emphasize  that  there  is  a  considerable 
redundancy  in  the  type  of  turnover  diagram  I  have  suggested;  i.e.,  it  contains 
less  information  than  is  apparent  at  a  first  glance.  This  is,  of  course, 
almost  necessary  because  otherwise,  from  what  I  just  said  in  the  previous 
paragraph,  there  would  be  no  limit  in  the  amount  of  "information”  one  could 
extract  from  bg(t)-data.  Part  of  this  redundancy  is  revealed  by  the  fact  that 
the  different  parameters  are  not  entirely  independent  of  each  other.  For 
instance,  it  is  easily  shown,  and  intuitively  obvious  that 

eu  -  *1-1 IV (ti-i>  -  W]  /  v°> 

i.e.,  the  size  of  the  initial  part  of  is  strongly  correlated  to  ^  By 

the  same  token,  0.o  is  positively  correlated  to  the  difference  It  -  t.  |  . 

1  i  l-l1 
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proved  in  the  first  reference  in  (2).  The  interpretation  in  terms  of 
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automatically  in  the  closed  logic,  which  the  physical  tracer  dynamics 
defines  (8). 

15.  The  time  of  sojourn  is  here  used  instead  of  the  more  appropriate  but 
clumsier  mean-time  of  sojourn.  It  should  not  be  confused  with  the  concept 
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